Phenomenological models of dielectric functions and screened Coulomb potential 
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The evolution of static dielectric constants and polarization with doping are analyzed and dis- 
cussed using established experimental results. The variation of screened Coulomb potential with 
doping is derived theoretically to justify the said evolutions. The latter justification arises due to 
the infiuence of ionization energy [Ex) on the static dielectric functions. The bare Coulomb po- 
tential and the Thomas-Fermi screened potential with finite carrier density were recovered when 
Ex ^ 00 and Ex ~* E% respectively. Basically, these phenomenological models are associated with 
the ionization energy based Fermi-Dirac statistics (iFDS). 
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1. INTRODUCTION 

Analogous to ferro- and paramagnetic systems that 
deal with spin alignments, Ferroelectrics (FE) ex- 
hibit spontaneous polarization below the Curie temper- 
ature {Tc) whereas the unpolarized state exists above 
Tc- The simplest Ferroelectric (FE) crystaL SnTe 
with NaCl structure was reported by Riedl et al. [ij back 
in 1965. The Ferroelectricity (FE) observed in SnTe 
is of displacive nature, identical with other Perovskite 
FE. Presently, the Perovskite structured (Ba,Sr)Ti04 
(BST) plays an enormous role for the development of 
technological applications due to its high spontaneous 
polarization, -IQ-^ Cm'^ at Tc « 290 K 0]. Al- 
though substituting Nb into the BST system reduces the 
twin formation and Tc significantly, its dielectric con- 
stant and spontaneous polarization was found to be in- 
creased sharply 0. Basically, FE are highly insulat- 
ing, which prohibits direct current transport measure- 
ments. Consequently, measurements on polarization and 
dielectric constant are of more relevant to characterize 
FE. Due to its dielectric effect, FE has been exploited 
in microwave devices such as tunable oscillators, phase 
shifters and varactors 0, |M IE 0- Moreover, FE also 
have the potential applications for detectors and de- 
vices based on the energy conversion principle namely, 
mechanical<^thermal<^electrical energies . 

Yet another technological advancement of FE is its 
functionality as a gate insulator in Field-Effect Transis- 
tor (FET) as reported extensively by de Boer et al. 
and Podzorov et al. ^3 using organic molecules namely, 
Tetracene, Pentacene and Rubrene. One of their interest- 
ing observation is the dielectric constant's inverse depen- 
dence on carrier mobility. Furthermore, developments 
of the dielectric gate that consists of La, Al and Oxy- 
gen to substitute the common SiO^ has been success- 
fully produced and even patented Perovskite FE 
with Ferromagnetic properties such as ( Tb,Y,Bi)M nO.'^ 
have been investigated experimentally |2, Il2l ll,*! Il4j due 



to its possible technological impact on magneto-electric 
media, including the Dynamic Random Access Memory 
(DRAM) devices. BST related oxides are studied in the 
thin film form, usually grown by pulsed-laser deposition 
technique. Parallel to this, there are reports on the in- 
fluence of the film's thickness with FE. In this respect, 
a generalized Ginzburg-Devonshire theory was employed 
to evaluate the behavior of the FE thin film coated with 
two metallic electrodes 15j. Zheng et al. 15] found that 
imperfect surface produces non-uniform polarization dis- 
tribution and subsequently gives rise to the phase transi- 
tion temperature or Tc. Apart from oxides, liquid crystal 
is one of the well known soft matter that works on the 
principle of molecular polarization |l6| . By applying an 
intense alternating electric field during Triglycine Sul- 
phate (TGS) crystal growth, Arunmozhi et al. |T3 con- 
cluded that this procedure leads to a longer relaxation 
time, which in turn reduces the dielectric constant and 
the spontaneous polarization. On the theoretical aspect, 
a thermodynamical model has been proposed to explain 
the T and E dependence of polarization and susceptibil- 
ity of TGS by Otolinska et al. 0|. Moreover, non-oxide 
alloy, TlInS2 was reported theoretically fl^ to explain 
the peculiarities in the dielectric susceptibility, which is 
connected with the assumption that there exists coexis- 
tence of proper and improper FE, somewhat mimicking 
the TbMnOa, which also demonstrates the said coexis- 
tence 0. 

Simple non-Fermi gas system and strongly corre- 
lated matter are known to exhibit remarkable elec- 
tronic properties with minuscule substitutional doping. 
Handling such systems theoretically with the ioniza- 
tion energy (fj) based Fermi-Dirac statistics (iFDS) 
has been shown to be precise. For example, simple 
systems such as the diluted magnetic semiconductors, 
Mno.02Geo.98 '2^1 and Gai-^Mn^rAs HHl^ as well as for 
complex crystals nam elv, Cuprate high-Tc superconduc- 
tors, YBa2Cu307 [ii Im I2I EiJiiJ^ and Manganite 
ferromagnets, Lai-ajCaj^MnOa l2lll23 have been shown 
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quantitatively to be within the scope of iFDS. There- 
fore, iFDS is further appUed here to scrutinize theoreti- 
cally the evolution of polarization, static dielectric con- 
stant and the screened Coulomb potential with doping in 
FE BST and its derivatives. The reasons to revivify the 
dielectric functions and the screened Coulomb potential 
with iFDS is two-fold. The first one is as stated above, 
while the second reason being, to further our evaluation 
on the actual role of phonon in the Cuprate supercon- 
ductors and also to justify why hole compensation effect 
is unnecessarily invoked in the First-Principles approach 
for the ferromagnetic semiconductors, both of which will 
not be addressed here. It is well known that for free- 
electrons, the Thomas-Fermi screening length determines 
the screened Coulomb potential while the bare Coulomb 
potential is recovered when the electrons density is liter- 
ally zero. However, it is highlighted here that the bare 
Coulomb potential can still be recovered in the presence 
of finite electrons density, provided that Sj approaches 
oo. In addition to the theoretical justifications, experi- 
mental results from Pauling jl^, Jaswal et al. [sOl and 
Bell et al. 'si'l will be used to further reinforce the accu- 
racy of iFDS. Unfortunately, the theory presented here 
is not yet suitable to predict the T dependence of dielec- 
tric constant and polarizability, rather the substitutional 
doping dependence of the mentioned parameters will be 
accentuated and analyzed in detail. 

2. THEORETICAL DETAILS 
2.1. Polarizability and optical dielectric function 

The electronic polarization, P of a particular crystal 
can be defined as 

P =C5]a,K-)E;-'. (1) 

3 

C, aj{nj) and E^"^*^' denote the atoms concentration, 
polarizability at atom j, which is proportional to elec- 
trons number in atom 7' and the local electric field at 
atom site j respectively [S^. The local electric field for 
individual ions can be written as Eiocai = E -f- jj^P by 
assuming a spherical geometry for individual ions and 
the crystals are purely ionic [23 ■ E is the macroscopic 
electric field while it— is the Lorentz factor. In order 
to understand the evolution of P with n as a result of 
doping, another equation that relates doping with n is 
required. To this end, iFDS is utilized that eventually 
gives 
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Here, /i = —XEp whereas Ej represents the electron- 
ion Coulomb attraction, which is also an accurate repre- 
sentation of doping-parameter ,21^ J2, 23^ .24^ 26^ . The 
derivation of iFDS, f{£i) = exp [ - - X{Einitiai state ± 
£1)] by employing the restrictive conditions, Ui = n 
and y^.T ( Ei n i tini stn.t^ ± £T)ini = E can be found in 
the Refs. 13, lillii, El • Etnuiai state dcuotcs the 
energy at certain initial state and the Lagrange mul- 
tipliers, He + X£i = -In [(n/F)(27rAftV™e)^^^] and 
fJ-h-XSx ^ In [(p/t/)(27rAn2/m,,)3/2] ^ y jg ^^le volume in 
k-space and Ep is the Fermi level at T = K. me.h is the 
electrons or holes mass, h = h/2Tr, h is the Planck con- 
stant, while n and p are the electrons and holes number. 
In the previous work [13,111, A = l/ksT since the varia- 
tion of n was with respect to the variation of T. Contrary 
to the previous objective, A in this case can be connected 
with the binding energy of an electron in a Hydrogen 
atom, which is given by E-n. — (TO/2fi.^)(e^/47reo)^(l/n^). 
Here, n denotes the principal quantum number. The rea- 
son is to calculate the variation of n due to electrostatic 
potential, of which the n as given in Eq. ^ is actually the 
induced electrons density, riind arises as a result of both 
external and induced electrostatic potential. To obtain 
such a result, one has to make use of the two restric- 
tive conditions introduced in iFDS as stated above that 
respectively give 

00 
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i.= i:^e---^ /k^expf-A^Vj^- 
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Solving the above integrals, one can arrive at i? = 
3n/2A = En., hence A = (127reo/e2)n2rB, after taking 
n — 1. rB represents the Bohr radius. Notice the es- 
sential factor, e^^^i''^'^ in Eq. Q that will be used to- 
gether with A to define the carrier density above Ep as 
a result of the external and induced electrostatic field 
and finally to explicitly derive the static dielectric func- 
tion and the screened Coulomb potential for the strongly 
correlated (non-Fermi gas) systems. The dielectric func- 
tion (efg;, k)) for an isotropic polarization can be defined 
as 13^ 
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e(c^,0) = l + ^5I"^K)Ef"'- (3) 
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e(w — 0) and e{uj — oo) are the static and optical 
dielectric constants respectively, while lo is the frequency. 
Again, one can easily notice the proportionality, e(0) cx 
aj{nj) oc e^*-^^^^^' without involving much triviality. In 
fact, the parameters, ajijij) and e(0) will be analyzed 
with respect to Sj at constant T and its correctness will 
be justified with the static dielectric function and also 
with experimental results. 



Since the standard linear approximation gives L(no) ~ 
Ep and dEp/dno = 2Ep/3no, one can rewrite Eq. © as 
given below 

dF° 

^ [n{x) - no] = Ef{x) - E% ^ e^{x). 
duQ 

n{x) -nn^ ^^eLp(x). (7) 

Utilizing iFDS or Eq. Q specifically, such that n = 
nind — n(x) — no, then Eq. |(7J) can be expressed as 



2.2. Static dielectric function 

If the positive ions in the background are allowed to 
have a sinusoidal variation that leads to the positive 
charge density, p'^{x) = nje + pea;t(k) sin(kx) in the 
X direction, then the response of the electrons charge 
density is p~{x) = —n'^e + pi„fi(k) sin(ka:). The term, 
Pextiy) sin(ka;) defines the external electrostatic field that 
acts on the electrons. In other words, the electrons will be 
deformed as a result of (/3ea;t(k) and Lpindiy)- The latter 
being the induced electrostatic potential. One can em- 
ploy the Poisson equation, V^(^ = —j^i '/'(k) = ipext(y) + 
^indi^) and p(k) = Pert(k) + Pmd(k) in order to arrive 
at m 



P(k) 
k eo 



(4) 



The FE here concerns with the displacement of charge 
(D(k)), thus e(0,k) can be defined HI as D(k) = 
e(0,k)E(k). Using, V • D = V ■ eE = pe:.t/eo, V-E = p/eo 
and the Poisson equation, one can also show that |32| 
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Now, at absolute zero (0 K), the Fermi level is given 
by Ep — 2^(37r^no)^^'^. However, in the presence of 

ifextix) and ipindix), Ep{x) = ■:^^[ZT:'^n{x)Y/^ . Us- 
ing the the Thomas-Fermi (TF) approximation, one can 
write Ef{x) - E% = e[(pext{x) + (pmdix)] = eip{x). As a 
consequence. Linearization can be carried out to estimate 
the induced carrier density, n{x) — no at certain point x, 
in which Ep(no) can be linearized as L(no) at no = n{x) 
to give 



dE° 

L(no) = Ef{x) + [no - n{x)\ 
dno ^ ■' 



(6) 



n = n{x) - no = ^rw^^^i^) [K^f - ^i)] ■ (8) 

Consequently, using limgj-^^oo "-(2;) = no, one can 
surmise that it is rather impossible to induce a re- 
gion of enhanced electron concentration. However, 
linig^^^o n(a;) = no[l -l- {3/2Ep)e(p{x)\ supports the 
free-electron or Boltzmann-particle systems, which also 
satisfies the original TF approximation. It is worth not- 
ing that, for the latter limit, the electrons with its total 
energy, TE > Ep are literally free, forming the Fermi gas. 
In contrast, Eq. © denies such scenario, as long as Ep < 
TE < £j. To see this effect clearly, the screened Coulomb 
potential is derived and discussed in the subsequent sec- 
tion. Employing Eq. Q and Pind{x) = ~[n^{x) — n^^e 
(as introduced earlier), the Fourier components of Eq. ^ 
can be written as 
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= -^e'o^p [ME'p-Sj)] J ^(k)e'k-dk 

— 00 
00 

■exp[A(i?° J P^e'^-'dk. (9) 



3noe^ 
'eo2E" 



Using Eq. Q with term by term division, the Fourier 
component of the static dielectric function at zero fre- 
quency, e(0, k) can be derived from Eq. © as given be- 
low 



e(0, k) = 1 + -f exp [X{E'^p - £x)] . 



(10) 



Firstly, notice that both e(0, k) and e(a;, 0) are propor- 
tional to the factor of exp [A(£'^ — £x)] and e(0, k) 
e{uj,0). Secondly, the TF screening parameter remains 
the same, /C^ = 3no/2eo£'^e^. Considering the Fermi- 
gas system with Ep — £1, then one can arrive at the 



original TF dielectric function, e(0, k) = 1 
however, £1 —>■ 00, then e(0,k) — > constant. 
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2.3. Screened Coulomb potential 

Imagine a test charge e is placed at certain en- 
ergy level (> Ep) in the background of a si- 
nusoidal variation of the positive ions, then the 
Fourier transformed electrostatic potential for the un- 
screened Coulomb potential, '/?<;(k) can be shown to 
be (j/k^eo using Eq. The Fourier inversion of 

^,(k) is ^,ir) = [l/{27r)^] 27Tk\q/eok')e^-^dk = 

[q/An^eo] d\<iJ^^e'^'-''°^U{cos9) = g/^neor as it 
should be. Subsequently, one can derive the screened 
Coulomb potential, (p{r) using Eq. ij^l, in which e(0, k) — 
(p,(k)/(p(k). As such, ^(k) = (g/eok2)[kV(kV/C2 ,)]. 
Notice the parameterization, IC^ j = K.^ exp [X{Ep —Sj)] 
that has been used for simplicity. Finally, the Fourier in- 
version of iy9(k) is (p{r) that can be derived as 



(p{r 



{2itY J e,{k^+Kl,) 



X / e*^''™"^d(cos6l) 



dk 



kdk 



27r^eorJ k^ + /C^ ^ 



■ sin(kr) 



exp 



(11) 



7rk denotes the area in k-space. The screened 
Coulomb potential in this case is a function of doping 
parameter, Sj and one can enumerate the evolution of 
f(r) with doping for non-Fermi gas systems. The £i 
that originates from iFDS can be used to justify that an 
electron to occupy a higher state N from the initial state 
M is more probable than from the initial state L if the 
condition £i{M) < £i{L) at certain T is satisfied. As for 
a hole to occupy a lower state M from the initial state N 
is more probable than to occupy the state L if the same 
condition above is satisfied. SmiUai state is the energy 
of a particle in a given system at a certain initial state 
and ranges from -|-cxd to for the electrons and to — oo 
for the holes. Simply put, the magnitude of f{r)M'^N 
due to the polarization initiated excitation of an elec- 
tron (hole) from the initial state M (N) to N (M) is 
quantitatively smaller than ip{r)L^N. This effect can be 
non-trivially verified from the Eq. . On the contrary 
however, placing a test charge with its TE > Ep and 
Ep = £j imply Fermi gas that eventually gives rise to the 



TF screening potential [33 given by 
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Add to 



that, Eq. also reduces to the unscreened Coulomb 
potential even at finite uq provided that £i oo. The 
limit, 5/ — > cxD conceptually means that the electrons in 
the presence of the sinusoidal potential of positive ions 



are infinitely rigid with zilch polarization and eventually 
can be thought of as a sphere of a net charge q, as seen 
by the previously placed test charge e. Therefore, the 
potential is indeed a bare Coulomb potential, q/Aneor. 
Refer to Fig. ^ to clearly observe the transition of the 
TF screened potential (•) to the Coulomb bare potential 
(■) as a result of iFDS, at finite carrier density (no ^ 0). 
The solid lines, • (£/ Ep) and ■ (£"/ oo) are ob- 
tained from Eq. Hll|l with its appropriate limits. 



3. DISCUSSION 

Having derived all the required functions, now it is 
possible to venture into the experimental data reported 
by various researchers starting from the 1920s. Notice 
that the above derivations strictly requires substitutional 
doping that significantly maintains a single phase or solid 
solution. Obviously, doping with a purpose of creating 
a second phase or multi-phases are not applicable with 
iFDS based theory. The Clausius-Mossotti relation in 
the optical range is given by [s^ 



Af^ + 2 



(12) 



Af represents the refractive index. Using this relation, 
the electronic polarizability of atom j, which is a func- 
tion of electron number in atom j (aj{nj)) for a wide 
variety of ions [2^ Is^l were calculated and tabled in the 
Ref. M- The magnitude of J^j c^jX^-j) can be assumed 
to be accurate since those ions (plasma) easily satisfy 
the Lorentz factor [2^ Is^ for spherical ions or atoms. 
Figure [21 a)-d) depict the relation between aj {rij ) and 
£1 for 1-|- ^ A+ ions respectively. Those experimen- 
tal data points [i^, HO] were fitted in Fig. [3L)-d) with 
ctjin-j) = X e:xjp[y], which is in accordance with the 
principle of Eqs. © and ©. Here, the fitting pa- 
rameters, X = constant of proportionality while y = 
X{Ep — £1). Both values {X and 3^), obtained from the 
fittings are given in the Fig. [2K)-d) itself. In accordance 
with iFDS, polarizing a 4-1- ion, say Ce^+ needs to over- 
come an energy proportional to 6325 kJmol"^ (the 5*'* 
ionization energy), while La^+, Ba^+ and Cs+ need to 
overcome the energies proportional to 4819 kJmol""'^ (the 
4*'' ionization energy), 3600 kJmol"^ (the 3'"'' ionization 
energy) and 2234.3 kJmol"^ (the 2"'^ ionization energy) 
respectively. The absolute values are actually equal to 
the energy needed to ionize an atom or ion such that the 
electron is excited to a distance r. However, the ioniza- 
tion energies stated above are for taking that particular 
electron to r — > 00. Considering this scenario, one can 
surmise that both X and y are predicted to reduce with 
increasing valence state. As anticipated, the decreas- 
ing magnitudes of both X and y with increasing valence 
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states have been calculated so as to fit the experimental 
data points. Furthermore, one can also understand the 
exponential decrease of aj{nj) with 1+ (Fig. 2+ 
(Fig. 12))), 3+ (Fig. 13;)) and 4+ (Fig. 01)) ions, which 
are as a result of increased Sj. Simply put, the expo- 
nential reduction of cxj{nj) for Cs+ ... Li+, Ba^+ 
... Be2+, La3+ -4 ... and Ce'^+ ^ ... C'^+ are due 
to increased £i. Similarly, the reduction of both X and 
3^ is also due to reduced aj{nj) for high valence state 
ions (4+) as compared with 1+ ions. It is important to 
realize that 3^ oc A, which implies that y increases with 
increasing ionic size since A cx n^r^ as defined earlier. 

Out of naivety, one should not assume that y is sup- 
posed to increase from Li"*" Be^+ C*"*" since 
the ionic size is increasing. This is because, different 
elements have different ionic properties due to different 
number of positive charge protons, which determines the 
Coulomb interactions. Actually, y will increase from {C, 
Si, Ti, Zr, Ce}4+ ^ {C, Si, Ti, Zr, Ce}^+ and so 
on for the respective ions. The average £i has been de- 
termined with, = J2i=i %" where z is the va- 
lence state. Prior to averaging, the 1"*, 2"'', 3'''^ and 4*'* 
ionization energies for all the elements mentioned above 
were taken from Ref. [s^l- In summary, the proportion- 
ality, aj{nj) oc e-^''^''~^'^ is indeed valid. On the other 
hand, the static dielectric constant, e(0, k) must also sat- 
isfy Eq. Iini). To this end, BaTiOg (•), SrTiOg (■), 
CaTiOa (a) and Bao.sSro.sTiOs (♦) samples obtained 
from Ref. |3J are plotted in Fig.|51as a function of average 
£i at T\ = Tc + 200 K. The average f/s for Ba^+ , Sr2+ 
and Ca^+ are given by 734, 807, and 867 kJmol~^ respec- 
tively. After applying iFDS, one can predict that e(0,k) 
should decrease from BaTiOs SrTiOs CaTiOa 
while the magnitude of e(0, k) for Bao.sSro.sTiOa is ex- 
pected to be between BaTiOa and SrTiOs. All these 
predictions are remarkably in accordance with Eq. (|10|l . 
In fact, those experimental data [3ll | have been repro- 
duced with significant accuracy using Eq. Hl()|l as indi- 
cated with a solid line. Equation Hl()|l has been rewrit- 
ten as e(0,k) = A'exp[y] in which the fitted values for 
X and y are given in Fig. |21 itself. Equation ((Tn|l has 
been rewritten in the stated form because e(0, k) 1 or 

^ exp \\{E^~£i)\ ^ 1 for BST ferroelectrics. As a con- 
sequence, one can write that X ~ K,1/\<i' and y remains 
the same as \{E'p — In other words, one can under- 
stand why the static dielectric constant becomes smaller 
exponentially from BaTiOs — > Bao.sSrg.sTiOa SrTiOs 
CaTiOa, which is due to increasing Ej (£'/(Ba^+) < 



4. CONCLUSIONS 

In conclusion, the evolution of polarizability, polariza- 
tion, optical-static dielectric constants and the screened 



FIG. 1: The screened electrostatic potential with iFDS effect 
has been calculated against r with Eq. pi^ . The solid lines 
are from Eq. lllll . In the limits, £i Ep and £i — > oo, one 
can recover the TF screened potential (l/re~'^'^) and bare 
Coulomb potential (l/r) respectively. Notice that l/re~'^''^ 
and 1/r are represented with filled circle and filled square 
respectively. 



FIG. 2: Polarizability as a function of electrons number, 
aj{nj) are calculated (solid lines) against £i in order to fit 
the experimental data points for a) 1+, b) 2+, c) 3+ and d) 
4+ ions using aj{nj) — Xexp[y], which is in accordance 
with the principle of Eq. ^ . The fitting parameters, X and 
y are a constant of proportionality and A(-E^ — £i) respec- 
tively. All the experimental data points given in a),b),c) and 
d) were obtained from Pauling '23| and Jaswal-Sharma [s^- 



Coulomb potential has been derived for non-Fermi gas 
system using iFDS. These functions are able to ex- 
plain the doping effect as well as reproduces the ex- 
perimental doping trend with high accuracy in the well 
known BST and its related Perovskites namely, BaTiOa, 
Bao.sSro.sTiOs, SrTiOg and CaTiOg. The relationship 
between polarizability and iFDS is also found to be highly 
precise as compared with Pauling's experimental data on 
1+, 2+, 3+ and 4-1- ions, which in turn justifies the ap- 
plicability of iFDS. 
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